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1. Introduction 



The Standard Model (SM) is very successful and has been tested by experiments with high 
precision. However, there still exist puzzles in SM from theoretical perspectives. Two famous 
ones of them are the family problem and the flavor hierarchy problem. In SM, 3 families of 
quarks and leptons have similar gauge interactions and it seems that two heavier generations 
replicate characters of the lightest generation; while the masses of quarks and charged leptons 
have the obvious hierarchical structure. The hierarchical structure also exists in the Cabibbo- 
Kobayashi-Maskawa (CKM) mixing matrix. SM accommodates 3 family fermions and their 
phenomenologies by adjusting the Yukawa couplings. It does not supply an interpretation for 
the origin of 3 families and the hierarchy structure of flavor parameters. 

There have been several different approaches to address these two problems. A natural 
and popular way is family symmetry. Froggatt and Nielsen [Q] suggested a horizontal U{1) 
symmetry to understand the hierarchical fermion mass structure; while continuous non-Abelian 
family symmetry, like SU{3), can also lead to very promising results both in the quark sector 
and in the lepton sector j^. Recently, triggered by the data from neutrino experiments, discrete 
non-Abelian family symmetries have attracted many attentions. It is found that the A4 family 
symmetry Q can produce the famous tri-bimaximal lepton mixing matrix Q , which is a good 
approximation to the best fit value of the neutrino experimental data 

Different from the family symmetry approach^ , the family and flavor problems are addressed 
from new approaches in extra dimensional framework. Due to the exponential behavior of 
fermion profiles in extra dimensions 0, the mass hierarchy of fermions can be produced by 
parameters of the same order naturally. While family problem also gets new interpretations in 
extra dimensions. Several groups of authors show that 3 families of SM in 4D can originate 



from 1 family in 6D 11 1 by making an appropriate gauge background or choosing a metric 



of special structure. Although there have been many progresses in these directions, there is still 
no a well celebrated model to address the family and flavor problems in the community. 

In this paper, we focus on the warped extra dimension approach to address the family prob- 
lem and the flavor problem. In the Randall-Sundrum warped 5D spacetime |12|, let fermions to 
propagate in the bulk. The profiles of fermion zero modes can be of the exponential behavior, 
which depend on the 5D bulk mass parameters. Due to this character, the fermion mass hier- 
archy can be reproduced by the 5D bulk mass parameters of the same order. This approach 
supplies a beautiful geometrical interpretation to the flavor hierarchy problem. In the paper 
P^ , we attempted to understand the origin of the same order 5D bulk mass parameters. The 
origin of the 5D bulk mass parameters has close relation to the family problem, because one 
mass parameter can stand for one family. We suggested a two-layer warped 6D spacetime, and 
begin with 1 family fermion in the 6D bulk. We reduce the spacetime from 6D to 5D at the 
first step. As a result, the 5D mass parameters emerge as eigenvalues of a ID Schrodinger-like 
equation (or Kluza-Klein (KK) modes in 5D). The spacetime metric can be chosen such that 
only 3 eigenvalues are permitted. So in this setup, 3 families in 5D can originate from 1 family 
in 6D. When we further reduce the 5D spacetime to the physical 4D at the second step, the 
zero modes of 3 family fermions in 5D produce 3 family fermions in 4D. By coupling with Higgs 



""^ However, for symmetry approach in extra dimensional framework, see 
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field in 4D, these zero modes get masses and produce tlie 3 family fermions in SM. According 
to these considerations, we can imagine that both the family problem and the flavor hierarchy 
problem can be addressed in such an approach. 1 family fermion in 6D produces 3 families 
in 5D. The 5D mass parameters are of the same order. The same order 5D mass parameters 
further produce hierarchical structure of fermions in 4D due to the exponential behavior of zero 
mode profiles. In the following part of this paper, we construct a specific model along the above 
ideas. We find that the numerical results of this model can be very close to the experimental 
data by adjusting parameters in this model. The mass hierarchies of quarks and charged lep- 
tons are produced. Supposing the neutrinos to be the Dirac ones, the small neutrino masses 
are obtained. The quark mixing matrix and the lepton mixing matrix are also very close to the 
realistic CKM matrix and Petrov-Maki-Nakagawa-Sakata (PMNS) matrix respectively. 

This paper is organized as follows. In section ||, we simply introduce the model suggested in 



1 13], and develop it to be a more realistic one. We notice several treatments which are different 
from that in our previous paper in section |2[ In section ^, we construct a model along 
the considerations in the introduction, and adjust the parameters of this model to show that 
the numerical results can be very close to the experimental data both in the quark sector and 
in the lepton sector. In section ^, we give some analytical treatments about our model. By 
these analytical treatments, we show that there are some common concise structures shared by 
quarks and leptons. We also make some further discussions about this model in section We 
make summaries in section ^. Several appendices are added for a clear understanding of the 
paper. 

2. Preliminaries for model building 

In this section, we introduce some necessary tools for the future model building in section ^ 



In subsection 2A, following the discussions in our previous paper |13], we introduce the basic 
setup to show that how several fermion families in 5D can originate from 1 family in 6D by 
KK decomposition in two-layer warped extra dimensions. The critical point is how we can 
get finite KK modes while there are infinite KK modes in the usual KK decomposition. In 



subsection 2.2, we introduce an example in the popular quantum mechanics textbook to show 
that finite bound KK modes can be obtained. This simple example helps us to understand 
the problem more clearly. In subsection |2.3| , we make more discussions about the characters 
of massive KK modes, and analyze an example which will be used in the model building in 
section 0. In subsection |2^, we reduce the 5D action to the 4D one, and display the necessary 



results for model building. In subsection we discuss a problem which makes the basic setup 



in subsection 2.1 to be not realistic. We further suggest a new setup to bypass the problem. 



This new setup will be used in the model building of section H, instead of the basic setup in 



subsection 2.1. In addition, we notice several treatments about the setup which differ from that 
in our previous paper |jl3|. We also notice these different treatments in several footnotes. The 
model building in section |^ largely depends on the discussions in this section. 

2.1 Introduction to the basic setup 

Following the discussions in [13|, we consider a 6D spacetime metric with the special two- layer 
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warped structure, 



ds'^ = B{zf [A{yf {'q^^^dx^'dx" + dy"^) + dz'^] . 



(2.1) 



We choose q^i, = diag(— 1, 1, 1, 1), and suppose that the two extra dimensions are both intervals. 
A massive Dirac fermion in this spacetime has the action 



S= d^xdydz^ \ ^ [^c^^F^Vm^ - Vu^e^r"^] - i m^^ \ , 



(2.2) 



where is the sechsbien, and Vm = Qm + ^^M^ab,^ab = zi^aj'^b] is the covariant derivative 
of spinor in curved spacetime. o and M = 0,1, 2, 3, 5, 6 stand for the flat indices and the curved 
indices in the tetrads respectively, m is a real number to ensure that the action Eq. ( p. 2D is 
hermitian. The gamma matrix representations are as follows, 
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, i = 1,2,3, 



(2.3) 



where /i = 0, 1, 2, 3 and are the Pauli matrices. 

As a first step, we reduce the 6D action Eq. ( |2.2| ) to the 5D one by KK decomposi- 
tions. Rewrite the 6D 8-component spinor with two 4-component spinors and expand these 
4-component spinors with 5D fields as follows. 



Xi 

X2 



/ n n 



Suppose that Fn{z) and Gn{z) subjugated to the conditions, 

^ + \b-^bA Fn{z) - mBFniz) + A„G„(z) = 0, 



d 
dz 



+ ^B-'B,^ Gn{z) + mBGn{z) - A„F„(z) = 0, 



(2.5) 
(2.6) 



where B^ 



^^j^. With the above conditions, the 6D action Eq. ( ^.2]) is reduced to 



S 



d xdy K„ 



Mr, 



d'^xdy MmnA^^itpm'(pn, 



dzB^ [F^Fn + G'^Gn 



dz {F^Fn + G'^Gn) 



dzB- 



rp* TP I 







/ 2 





A™ + Xfi 



where we have defined the transformations 

Fn{z) = B{zY''I^Fn{z), Gn{z) = B{z)-^'^ Gn{z) . 



(2.7) 
(2.8) 
, (2.9) 

(2.10) 
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The above equations are satisfied for all KK modes, including zero mode and massive modes. 
By the transformations Eq. ( |2.10| ), equations ( |2.5| ) and can be simplified to 

^ - mB^ Fn{z) + XnGniz) = 0, (2.11) 

^ + mB^Gn{z)-XnFn{z) = 0. (2.12) 
Our purpose is to obtain the conventional effective 5D action 

- / d^xdyA\XnMn- (2.13) 



We consider two cases: 

Case (I): The different KK modes are orthogonal, that is, the normalization conditions 

Kmn = j dz (F4F„ + G*^Gn) = (2.14) 

are satisfied. In these conditions, the 6D action Eq. ( |2.7D reduce to be the 5D action Eq. ( |2.13| ) 
naturally. 

Case (II): The second case is that the normalization conditions Eq. ( p. 14 ) are not satisfied. 



In this case, K and M are both matrices. It seems that we can not obtain the conventional 
effective 5-dimensional action Eq. ( |2.13D at first sight. However, if the number of KK modes 
is finite and the matrix K is positive-definite, we can redefine the fermion fields to obtain an 



action, which has the same form with that of Eq. ( 2.13 ). The difference is that the eigenvalues 
An are modified to different values. The redefinition process can be done in two steps. 
Step (1): Decompose the hermitian matrix K as follows, 

K = yUy = h'^h, h = y/Xv, (2.15) 

A = diag(Ai,A2,-- - ,A„), 
VA = diag(v/A7, yX^, • • • , y/K). 

In the above expressions, Aj > 0, z = 1, 2, • • • , n, as we have supposed that K is positive-definite. 
Redefine ■0n as 

Tpm = Hmni'n, (2.16) 

then in the new basis ^l^n, M becomes 

M = {H-'^)''MH-\ (2.17) 

After this step, the kinetic term of Eq. ( p.7| ) becomes the conventional form as that in Eq. (|2.13| ); 
Step (2): Obviously M is a hermitian matrix. We can diagonalize this new matrix as 

M = U'^AU, (2.18) 
A = diag(Ai, A2, • • • , A^). 
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Redefine also the new basis 'ipn, 

ipm = Umn'4'n- (2.19) 

We can obtain the action 



-Y. I d'^xdyAhXni^Jn, (2.20) 



which has the same form with the action Eq. ( |2.13| ). The total operations on the fermion fields 
equal to 

^ =Uij = UH'4j, (2.21) 
= H~^^ = H~^U^, (2.22) 

where we have omitted the subscripts. 

In the above, beginning with the 6D action Eq. (|2.2| ), we obtain the 5D action Eq. ( 2.13| ) 
by KK decompositions. We can interpret this process as follows: the action Eq. ( p.2[ ) stands for 
1 fermion family in 6D; while the action Eq. ( 2.13| ) or Eq. ( |2.2[1| ) can stand for several fermion 
families in 5D if we can restrict the number of the KK modes to be finite. The masses of these 
KK modes are determined by the equations ( p. 11 ) and ( |2.12| ). Moreover, it is not necessary to 
require that these masses have the hierarchical structure. These 5D bulk masses of the same 
order are enough to produce the hierarchical structure of 4D fermions according to the works 
in P|. Therefore, the key point is how we can get finite KK modes. This issue is the theme of 
the following subsection. 

2.2 An example for finite bound KK states from ID quantum mechanics 

Now we analyze the solutions of equations ( 2.11| ) and ( |2.12| ). For a zero mode (A = 0), these 
equations decouple and are easy to be solved. The solutions are given by 

Fo{z) = ^^exp ( [ mB{C)dc] or 0,Go(^) = ^^exp (- [ mB{C)dc] or 0. (2.23) 

Vf^Vo \Jzo J JiAT \ Jzo J 



We introduce I of the length dimension in order to make the normalization constants to be 
dimensionless. For massive modes, we can combine the first order differential equations to 
obtain second order equations 

-^Fn{z) + [-mB, - m^B^] Fn{z) + \lFn{z) = 0, (2.24) 

^Gn{z) + [mB, - m^B^] G„(z) + X^Gniz) = 0. (2.25) 
Rewriting them in another form, we see that they are similar to the ID Schrodinger equations 

- ^Pn{z) + V{z)Fn{z) = XlFniz), (2.26) 

- ^Gn{z) + V{z)Gn{z) = XlGniz), (2.27) 
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with potentials 



V{z) 



mB^+m^B'^, V{z) 



-mBz + m^S^, 



(2.28) 



where B^ = — ^j-^- In order to obtain the effective 5D action Eq. ( |2.13| ), we should require that 
the integrands in Eq. ( |2.8D and Eq. ( p. 9] ) are finite. These requirements can be satisfied if the 
KK modes are the bound states of Schrodinger equations ( 2.26| ) and ( ^.27 ). So the problem 
how we can get finite KK modes transforms to the problem how we can get finite bound states 
of the ID Schrodinger equations ( |2.26|) and (|2.27D . 

There is a simple example |14| in ID quantum mechanics which has finite bound states. It 
is the square potential well of finite depth and width. The potential V{z) is given by 



V{z) 



.|>f, 



(2.29) 



where V > The Schrodinger equation is 



— -^u{z) + V{z)u{z)=Eu{z). 
2m az^ 



(2.30) 



The existence of bound states requires that E < V. We first consider the case E > 0. The 
solutions of Eq. ( 2.30 ) can be classified by parity. 
The solution of odd parity is given by 



u{z) 



—C exp(/3z) 
B sin(az), 
C exp(— /3z) 



z < —- 



Z \< - 



(2.31) 



z > 



where a 



2mE 



2m{V-E) 

— w — • 



The continuity of u{z) and ^^^^ at z = ib| requires that 



- ^ cot ^ 



V, 

mV 
2F' 



(2.32) 
(2.33) 



where ^ = a|, rj = (3^. 

For the solution of even parity, we obtain 



u{z) 



C exp(/?2;), 
B cos(a2;), 
C exp(— /3z) 



z < -- 



z\< - 



Z > 



2 ■ 



(2.34) 



The continuity conditions require that 



^tan^ 



V, 

mV 
2F' 



(2.35) 
(2.36) 
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The eigenvalues are determined by equations ( ^.32 ), ( 2.33| ), (|2.35 ) and ( |2.36| ). Their number is 
finite^, and the condition that there exist the very n eigenvalues is that 



From Eq. we know that the number of eigenvalues is restricted by the depth and the 



width of the square potential well. 

In the above, we have discussed the situation E > 0. However, in the process of KK 
decomposition, it is not necessary to require that E > 0. We should also discuss the situation 
E < 0. We can obtain solutions of < by replacing a and ^ with ia and in the solutions 
of > 0. We display only the conditions that determine the eigenvalues here. For the odd 
parity solution, we have 



-|coth| = r?, (2.38) 
mV 



~r) 9 mV 9 

+^ = ^« ' (2-39) 



and for the even parity solution, we have 

-^tanh| = r/, (2.40) 



where a = y ~^^) ^ = ci;f j while (3 and rj keep invariant. It is obvious there are no solutions 
that satisfy the above equations. 

In summary, there are finite bound states in the square potential well. This is a simple 
example to show that it is possible to obtain finite bound states by choosing the potential 
properly. We are also interested in what kind of metric B{z) can make the square potential 
well. Replacing V{z) with the square potential well in Eq. ( ^.28]) , we obtain 

{0 \ z\< - 

' ' 2' (2.42) 
V, Ul>f. 

Solving it for B{z), we obtain 

'' ' ' z\< - 

^ 1^ 9) 



^tanhiVVz + C), 



z \> - 



where C and C are constants. They can be determined by boundary conditions and continuity 
conditions. This metric is very similar to that in the papers jl^], where the authors opened up 
the extra dimension at infinity to address the cosmological acceleration problem. For another 
example for finite bound solutions in different background, see [llf 



For details, see 
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2.3 General characters of massive modes and an example for model building 

In the above subsection, we give an example in ID quantum mechanics to show that it is possible 
to obtain finite bound modes. In this subsection, we give another example which also has finite 
bound modes. This example will be used in the model building of section |3[ Before doing that, 
we analyze some general characters of equations ( ^.ll| ) and ( 2.12 ). 

The zero mode solution of equations ( 2.11 ) and ( p. 121) has been given in Eq. (112^) . The 
massive mode solutions are determined by the second order equations ( 2.26 ) and ( 2.27| ). From 
equations (|]2|) and (|23^ ), we see that eigenvalues A„ emerge with the form A^; so for =bA„, 
equations ( ^lej) and ( ^ ) are both satisfied. In fact, we can infer from equations ( p.ll[ ) and 
( |2T2|) that if the ])air (F„,G„) is a solution of equations (|2.11| ) and ( p.l2|) corresponding to 
the eigenvalue A„, then another pair (F„, — G„) is also a solution of equations ( 2.111) and ( 2.12| ) 
corresponding to the eigenvalue —An. Therefore, the massive solutions always emerge in pairs^. 
We denote these pair solutions explicitly as follows, 




~ An 



Fn 



-Gn 



(2.44) 



Therefore, if these massive modes are discrete, we obtain (2n + 1) modes in all, where n = 
0, 1, 2, 3, • • • , including the zero mode and massive modes. Supposing that there only exist a 
pair of massive modes, we can write the spectrum of equations ( 2.11 ) and ( 2.12| ) according to 
above discussions as 




spectrum of A — > { . (2-45) 

Now we analyze an example which will be used in the model building of section ^. As in 

j) to be 

s, a,b,uj > 0. 



our previous paper [13|, we choose the metric B{z) to be 



^(^) = «^^IJ7TT' s,a,b,u;>0. (2.46) 

We choose the extent of z to be the semi-infinite interval [R, oo). This metric is similar to 
that we obtained in Eq. ( |2.43| ). We implicitly suppose m > in following discussions, unless 
we announce it explicitly. 

For this metric, the normalizable zero mode solution is given by 



Fo{z) = 0, Goiz) 



Wo' VT 



(2.47) 



For massive modes, Eq. ( p.26| ) can be solved by hypergeometrical functions, 
F[z) = Cxe'^^^e"^^ + 6)^~''hypergeom ( p - /i + z^, 1 - p - + z^; 1 - 2/u 



+ C^e^'^^ef''' + 6)"^"''hypergeom p + /i + z^, l-p + /i + i/;l + 2^ 



gW2 + h 



(2.48) 



^In our previous paper we only keep the zero mode and the positive massive modes, and omit the negative 
massive modes. 
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where p = ^s{l - fi 



and V 



Ci and C2 are 



constants. We have omitted the subscript n explicitly. We only give the solution for F{z) 
here. The solution for G{z) can be determined by F[z) through Eq. ( 2.11| ) or by Eq. ( |2.27| ) 
directly. As discussed in jl^, in order to make the solution to be finite when z — > 00, the 
hypergeometrical series must be cut off to be a polynomial. This cut off can be completed by 
four different kinds of choices. We can cut off the series after the coefficient Ci by conditions, 



p- fj, + u = -n, 
or 1 — p — /x + i/ = — n, 



(2.49) 
(2.50) 



in which n = 0,1,2,3, 
conditions, 



Otherwise we can cut off the series after the coefficient C2 by 



p + p + u = -n, 
or 1 — p + p + i' = —n. 



(2.51) 
(2.52) 



The four equations (|2.49|) , (|2.50|) , (|2.51|) and ( p.52|) are not necessary to be satisfied at the same 
time. Anyone of them is enough to cut off the series. So we can obtain four kinds of solutions 
corresponding to these four different cut off ways generally^. The next important task is to 
analyze whether the number of eigenvalues determined by these equations are finite. We discuss 
this problem in appendix In this subsection, we only display the results. 

In this paper, we are interested in the situation f > 1. According to discussions in appendix 
we can have two different ways (2.49) and ( |2.50| ) to cut off the series. We might obtain two 

kinds of solutions corresponding to these two cut off ways. While equations ( p. 51 ) and ( 2.52| ) 

have no solutions, so the corresponding cut off ways do not work. As we analyzed in appendix 
n must be finite in equations ( 2.49| ) and ( 2.50| ). The extents of n are given by equations 

(|Al^ ) and ([All . For convenience of reading, we copy equations ( |A.13D and ( [A.14| ) with new 

sequence numbers here. 



m f a \ m 
— s { - — 1 < n < — s 



m / a \ m 
— s --l <n+l< —s 



2 /a 



for (12:491) , (2.53) 
for pISOl). (2.54) 



From (^^5^ ) and ( ^Isl ), we know that the number of eigenvalues is determined by a pair of 
parameters (|, ^s). We can choose the values of this pair of parameters to make the very two 
massive modes left. Note that conditions ( 2.53| ) and ( 2.54| ) are not necessary to be satisfied at 
the same time. They are different restrictions for two kinds of cut off ways respectively. 



*In our previous paper we only consider the solutions corresponding to the cut off condition Eq. (2.5C) 
The other three kinds of solutions are omitted. 
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2.5 




Figure 1: The parameter set for 3 families. 



The values of —s) construct a 2D plane. A point set in this 2D plane restricts n to be 



a specific value. We give a point set of 



ps) determined by following equations 
11 < 1, 



ni 



ijj 



m 
1 < —s 

to 



m 
— s 

UJ 



1 + 



1 



a\ 2 



a 



< 2, 



< 1. 



(2.55) 
(2.56) 

(2.57) 



Conditions ( 2.55| ) and ( 2.56| ) imply that only n = 1 is permitted in Eq. ( p. 53 ); so there is one 
solution corresponding to this kind of cut off. The condition (|2.57|) implies that Eq. (2.54) is 



impossible; so there is no solutions corresponding to this kind of cut off. As discussed in the 
above, the massive modes always emerge in pair. So this parameter set is enough to make a 
pair massive modes. Together with the zero mode, we obtain the very 3 modes in all. For the 
future model building, we give solutions for these modes explicitly in appendix |^. 

The point set determined by equations ( |2.55| ), ( |2.56| ) and (|2.57 ) forms a specific area in 2D 
plane. This area can be visualized in Figure. |l|. 

For X > this area is given by the area between the curve v = ... } , — — and the 



curve y 



y 



this area is given by the area between the curve 



and the curve y 



However, the points at the curve y = 



must be excluded. As we discussed in appendix]^, in this curve, the equation 



1 



m /a 

UJ \0 



1 



(2.58) 



is satisfied. According to Eq. (|A.2l| ), when Eq. (|2.58|) is satisfied, A = for Eq. (|2.49|) . So 
in this situation, n = 1 gives a massless solution but not a massive solution. It is not a new 
solution and it coincides with the zero mode solution given by Eq. ( ^.47 ). We should exclude 
the points at this curve in order to ensure that we obtain massive solutions. 
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2.4 Reducing spacetime from 5D to 4D 



We have reduced the spacetime form 6D to 5D in subsection 2.1. In this subsection, we further 
reduce the spacetime form 5D to 4D. We begin with the effective 5D action ( |2.2[1| ). Note that 
masses A„ are real irrespective of values of in equations ( 2.11 ) and ( |2.12| )^, because the 
induced matrix M in (2.17) is hermitian. As in the interval approach of RS model, we choose 
the extent of y to be a finite interval [L, L']. 

We make KK decompositions by expanding the 5D field V' with the 4D chiral fields as 

^ = Y,A{y)-\f'i\y)i;^l^ + f'i\y)^P%\ (2.59) 
i 

where ipL = —^^ipL, '^R = 7^V'-R- Iii the above expanding, we have omitted the subscript for 
ip. It is understood that we do the similar operation for each tpn in the action Eq. ( |2.20| ). The 
KK modes satisfy following equations 

+ A^) f'i\y) - Mi^^^fii\y) = 0, (2.60) 
a - A^) fii\y) + ^^'^f?iy) = 0, (2.61) 

where m^-'^ are eigenvalues. 

In this paper, following a popularly adopted approach in papers ^J, we employ the zero 
mode approximation (ZMA) approach, that is, we solve for the fermion bulk profiles without the 
brane interaction terms at first and then we treat the brane interaction terms as a perturbation. 
In this ZMA approach, we can choose two groups of consistent boundary conditions for equations 
and ([I^ ) as follows. 

Group (I) : + A^) /^(y) = 0, fl(\y) = 0, at y = L and L' , (2.62) 

Group (II) : (^^ - XA^ fi(\y) = 0, fi^\y) = 0, at y = L and L' . (2.63) 

The boundary conditions ( 2.62] ) save the left-handed zero mode and kill the right-handed one; 
while the boundary conditions ( 2.63| ) do the opposite. 

In the model building in section ^, we will choose the metric factor A(y) to be 

A{y) = -, (2.64) 

y 

that is, the RS spacetime. By solving the bulk equations for this metric, we obtain the normal- 
ized left-handed zero mode profile and the right-handed one as 



^In our previous paper we require that A,i must be real. It is a superfluous requirement and is not 
necessary. As in the discussions in appendix H, X„ can be pure imaginary. 
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where c = AL, e = jj, and we suppose that e <C 1. In this paper, we will also use the canonical 
normalized zero mode profiles and their values at y = L' . Their values at y = L' are given by 



f(0) 
JL 

m 

JR 



2c- 1 



c2c-l 



e 2, 



-2c -1 

-e 2. 



C-2C-1 



(2.67) 
(2.68) 



where fi\{y) = V^fi^iv)- For discussions in the next subsection, it is useful to find the 
asymptotic behavior of these profiles. 



VI - 2c, for c < i, 

k' for c ^ i, /i?^ ~ < 

V2c- 1 e'^"^, for c> i. 



-Ine' 



V-1 - 2c 



for 
for 
for 



c< i, 

c^i, (2.69) 



We have supposed that e <C 1 in above expressions. 

In the above, we discuss characters of zero modes. In the following we discuss massive 
modes. For massive modes, equations ( 2.60 ) and ( 2.6l| ) can be combined to be the second order 
equations 



dy^ 
a2 



dy- 



:fij\y) + [-A^, - X'A']fil\y) + :m(^^'f^^\y) = 0, 



(2.70) 
(2.71) 



where A. 



For the metric (2.64), equations ( |2.70| ) and ( p. 711 ) can be solved by Bessel 



dy 

functions. Following the standard procedure of analyzing Sturm-Liouville equations, we can 
obtain from equations ( |2.70| ) and ( |2.7l| ) that 



m 



(j)2 



m 



(i)2 



(m«2) 



(m»2)* 



dy(ftfi^^ 



Ai) d M)* 

JR j.Jb 



SR ^ySR 



d 



dy 



IR 



iJ) 



(2.72) 
(2.73) 



The above equations are derived from equations ( p.70| ) and (2.71) for massive modes, but it 
also applies for zero modes. The boundary condition ( 2.62| ) ensures that left-handed and right- 
handed KK modes satisfy orthogonal conditions 



(2.74) 



While the boundary condition (|2.63| ) also makes orthogonal conditions (|2.74| ) satisfied. 



2.5 A more realistic setup for model building 



According to discussions in subsection |2.3| subsection 2^, we find a problem which prevents 
us to construct a realistic model. This problem arises from the following contradiction. On 
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one side, from the asymptotic expressions of the zero mode profiles ( 2.6S| ) in subsection 2A 
we know that in order to produce fermion mass hierarchies, the zero mode profiles should be 



of exponential behavior. According to Eq. ( 2.69 ), because e <C 1, we need positive parameters 

for right-handed zero 



-c > 



c > ^ for left-handed zero mode profile, and negative parameters 
mode profile. On the other side, the parameter c comes from the 5D mass A. We have found 



in subsection 2.3 that the spectrum of A consists of one zero mode and two massive modes itA 
in pairs^. For left-handed zero modes, A = and —A prevent the profiles to have exponential 
behavior; while for right-handed zero modes, A = and A prevent them to have exponential 
behavior. This contradiction constructs an obstacle to produce the fermion mass hierarchies. 
In this subsection, we construct a new setup to bypass this problem. 
Instead of the action (|2.2| ), we suggest a new bulk action as follows, 

S= /d^xdydzV^ [^e^r'^VTv^- Vjv^e^r"^'] -i m^^'-A^^T^eljr^^'l , (2.75) 



where is defined by Eq. (|2.3| ). Our metric Ansatz (2J) keeps invariant. Ai is real to ensure 
that this new action is hermitian. Comparing it with the action (|2.2D , we add a new term. 
This new term has a strange form, as it breaks the 6D local lorentz invariance obviously. Here 
we regard this new action as an effective one to bypass the problem discussed above. We will 
discuss the possible origin of this new action in section ^. 

Making the KK decompositions as in subsection 2.1, equations (2.11) and ( p. 12 ) are mod- 
ified to be 



d_ 

dz 
d 
dz 



- mBj Fn{z) + (A„ - M)Gn{z) = 0, 
+ mB] Gn{z) - (A„ - M)Fn{z) = 0. 



(2.76) 
(2.77) 



The induced 5D action is the same with that in equations (^j) , (12. 8| ) and (|2.9| ) . The differences 
are that the KK modes and the eigenvalues A^ determined by the new equations ( |2.76| ) and 
( 2.77 ) instead of equations ( |2.11| ) and ( ^.12| ). By defining new variables 



we obtain 



AA„ = A„ - 7W, 



— -mB]Fn{z)+A\nGn{z) 



dz 



+ mB ] Gn{z) - AXnFniz) 



0. 



(2.78) 

(2.79) 
(2.80) 



These new equations have the same form with equations ( 2.11 ) and ( ^.12 ). So if we replace 
Xn with AA„ in equations ( 2.11| ) and ( ^.12 ), the results that we obtain in subsections 2^2^ , 



6 



In subsection 



2.3, we show that the spectrum of A consists of one zero mode and two massive modes ±A in 



pairs. However, we can prove that this resuhs also applies to A. We wiU give the proofs in subsection 4.1. Or 



one can employ the example in subsection 2.3 to check it with numerical method directly. 
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appendix ^ and appendix ^ apply here. Of course, we should notice that in Eq. (|2.78| ) 
emerge in actions (^), (2.8) and ( ^ ) but not the new variables AA„. We will give more 
details in section There is a zero mode solution corresponding to AA = 0, and the massive 
mode solutions also emerge in pairs. A pair of massive modes for AA^ are given by 



AA„ 




AA„ 



-Gn 



(2.81) 



The spectrum of AA will be the same with that of A given by Eq. ( p. 45 ). While by Eq. ( 2.78| ), 
the spectrum formula (|2.45|) of An now changes to be 



spectrum of A A 



AAi 

A„ = AAn + spectrum of A 
> 

-AAi 



AAi + 7W 

M .(2.82) 

-AAi 



By the new spectrum Eq. (|2.82[) , we can bypass the problem discussed at the onset of 
this subsection. From equations ( p.79 ) and ( 2.80| ), we know that AA„ depends on the mass 
parameter m and the metric B{z). So just like we discussed in subsection |2.3| , the number and 
the spectrum of AA^ are completely determined by m and parameters in B{z). The number 
and the spectrum of A A are irrelevant to the new parameter M. We can adjust M to change 
the spectrum of the final eigenvalues A. If the value of 7W is a large positive number, then the 
spectrum of A can be all positive by Eq. (2.82), which are expected for left-handed zero modes. 
We can also make the spectrum of A to be all negative by choosing small negative value of A4, 



which are expected for right-handed zero modes. As we will discuss in subsection |4.1| , the new 
spectrum Eq. ( ^.82 ) for A also applies to A. The above discussions for A also apply to A. These 
features bypass the problem at the beginning of this subsection. 

In the model building of section ^ we will adopt this new setup in this subsection. Com- 
pared with the setup in subsection |2.1|, it is not more difficult to analyze this new setup. We 



only define the new variable by Eq. (2.78). All other analysis in subsection 2.2, subsection 2.3 
appendix ^ and appendix ^ apply here. 

3. A model of fermion mass and mixing 



In subsection 2.1, we discuss how several families in 5D can generate from 1 family in 6D. In 



subsection 2^, we suggest a new setup. By this new setup, we saw that the fermion mass hier- 
archies in 4D can be produced when we further reduce the spacetime from 5D to 4D. Therefore, 
we hope that fermion mass hierarchy problem and family problem can be both addressed in 
such an approach. In such an approach, the mass hierarchy problem will closely relate to the 
family problem. In this section, we construct a specific model along these discussions. In sub- 
section p.l|, we introduce the model in the 6D bulk spacetime. In order to construct a realistic 



scenario, we make two assumptions in the model building. In subsection |3.2| , in order to obtain 
a 4D effective theories, we employ two step KK decompositions to reduce the spacetime from 
6D to 4D. At the first step, we reduce the spacetime from 6D to 5D, then we further reduce the 
spacetime from 5D to 4D. After deriving the 4D effective theories, we further give numerical 
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results in subsection p.3|. In subsection 3.3.1, we give numerical results in quark sector. In sub- 



section |3.3.2| , supposing neutrinos to be Dirac ones, we apply this model to the lepton sector. 
In both sectors, the numerical results can be very close to the experimental data. We also give 
brief comments about the numerical results in subsection 3.3.3| . 



3.1 Model building 

We begin with 6D spacetime with the two-layer warped metric 

ds"^ = QMNdx^dx^ = B{zf [A{yf [T^^^dx^'dx" + dy"^) + dz^ 
We also choose the metric factors to be 

, My) = -■ 

y 



B{z) 



(3.1) 



(3.2) 



The parameters in this metric will be designated when we give numerical results in subsection 
3.3. As discussed in section ^, we choose the two extra dimensions are both intervals. The 
extent of z is a semi-infinite interval [i?, oo]; while the extent of y is a finite interval [L, L']. By 
these choices, the spacetime is sandwiched by 3 co-dimension 1 4-branes sited at z = R, y = L 
and y = L' . We also introduce two co-dimension 2 3-branes sited at the brane intersections 
{z = R, y = L) and {z = R, y = L'). They can be dubbed UV-brane and IR-brane respectively. 
We will only designate the field contents on IR-brane sited at {z = R, y = L'), while the field 
contents on UV-brane are omitted in the following discussions. The metric on these 3-branes 
are given by the induced metric 



g^,y{x,R,L), 
9l^{^) = 9^w{x,R,L'). 
We discuss the quark sector at first. We introduce the quark field contents as 



U 



V. 



(3.3) 
(3.4) 



(3.5) 



They transform under the gauge group SU(3)c x SU(2)l x U(l)y as Q = (3, 2)_|_i/6, U = 
(3, l)+2/3i ^ = (3, l)_i/3. These fields are 6D fields. Note that there are no family indices for 
them, that is, we introduce only 1 family in the 6D bulk. Their actions are given by 



Sr 



Su 



Sv 



d'^xdydzyj —g 



— / d'^xdydz\ 



I 



1 &yan 



d xdydz^J —g< - 



d'^xdydz^f^MuUV'^ elV'U 



d xdydzy/^\ - 



imuUU 



imj}T>T> 



d'^xdydz^/^MvW' elV'V. 



(3.6) 



(3.7) 



(3.8) 
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As in the usual field theory, we introduce the interactions between fermion field and gauge 
filed by requiring that actions ( |3.6| ), ( |3.7D and ( |3.^ ) are invariant under the 6D local gauge 
transformation 



U 



^-i[gej{x,y,z)S^+g'e{x,y,z)](^ 



(3.9) 



where ^ are the generators of the gauge group SU(2)j^''. However, we do not employ this 6D 
local gauge transformation ( |3.9| ) in this paper. Instead of it, we assume that the fermion actions 
are only invariant under the 4D local gauge transformation 



^-i[g§^(x)^+g'§(x)]q^ 
g-ig'e(x)^ p ^ ^-ig'§{x)jy^ 



(3.10) 



Here the gauge parameters only depend on the 4D coordinates. Requiring the invariance under 
transformations Eq. ( |3.10| ), we obtain the interaction terms 



•^Qint 
Suint 
Svint 



gW^,{x)-+g'B^{x) 



d'^xdydzV^ {lie^r''[g'B^{x)]U} , 
d^xdydz^ {Pe^r"[5'S^(x)]P} , 



(3.11) 
(3.12) 
(3.13) 



where fi = 0,1,2,3. Because we suppose only the 4D local gauge invariance, the gauge field 
components W^, Wq and B^, Bq relevant to extra dimensions are not necessary to be introduced 
here. These gauge fields only depend on the 4D coordinates. We designate their actions to be 
the brane actions on the IR-brane 



•Sgauge = / d'^xdydz^/^g^ 



9 vis 9 vis 



1 ^ 

4' 



5{y-L')5{z-R)\, (3.14) 



where g^^^ is given by Eq. ( |3.4D . 

Note that interaction terms ( |3.11| ), (|3.12 ) and ( |3.13| ) are different from conventional brane 
interaction terms 



Si 



Qint 



d^xdydz,/^,{UeZra[g'B^{x)]l(5iy - L')5{z - R)] 



(3.15) 



where e"(l, a, ;U = 0, 1,2,3 are the tetrad determined by the IR-brane metric g^'^ . Interaction 



terms ( 3.111 ), ( p. 12 ) and ( p. 13 ) are introduced as above to ensure that we can obtain an unitary 
mixing matrix in the ZMA approach; while the brane interaction terms like Eq. ( p. 15]) break the 
unitarity of mixing matrix remarkably in our present model. Interaction terms ( |3.1l| ), ( |3.12| ) 
and (|1|) are our important assumptions in this paper. We can also introduce gauge fields in 



^In this paper, we focus on the electroweak interaction sector, so the color indices are imphcit. We can 
introduce the color interaction just like that we do for electroweak interaction. 
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the 5D bulk as in the papers [0] . The unitarity of mixing matrix in the ZMA approach is also 
kept. However, in this paper, we employ the more simple approach as introduced in the above. 

Now we introduce the Yukawa interactions between fermion fields and the Higgs fields. We 
designate them as 

^T^Yukawa = J d^xdydz^ S^^iQV + /Jj^QF^P] $(x) - L') + H.C.j , (3.16) 

^WYukawa = j d'^xdydzV^ l^^iQU + f3uQT'UMx)^5iy - L') + H.C.j , (3.17) 

V2 \ h[.x) + i(t)i[x) J V2 y-0i(x) +z02(a;) J 

where $(x) is a doublet of SU(2)i. A is a constant with mass dimension. The factor ^ emerges 
because our metric Ansatz Eq. ( |3.1j ) has the conformal form. If we employ the Gauss normal 
coordinates for the sub 5D metric, the factor is superfluous. The parameters fu and fv 
are real numbers and dimensionless; while I3u and /3x> can be complex numbers generally. We 
also introduce the terms after f3ii and Z?©, which break the parity symmetry of the action. As 
we will discuss in following subsections, these terms are the origin of CP violation. If we drop 
these terms, the effective 4D actions will be CP invariant measured by the Jarlskog invariant 
measure. 

Note that the Yukawa interactions are different from the conventional brane Yukawa inter- 
actions 

SvBr.ne = j d^xdydz^^ | ^ [QP + /J^^QT^P] $(x) ^% - L')5(z " i?) | 

+ H.C., (3.19) 

^WBrane = j d^xdydz^/^ | ^ [QZ^ + (3uQT^U]^x)^5{y - L')5{z - i?)| 

+ H.C.. (3.20) 

We adopt interactions ( |3.16| ) and ( |3.17] ) instead of the brane interactions ( |3.19| ) and ( |3.2[1| ). 



The reasons are as follows. By the numerical method in section 3.3, we found that interactions 



(3.19) and ( 3.2C| ) lead to a massless fermion in the ZMA approach in the 4D effective theories. 



While interactions ( |3.16| ) and ( ^.17 ) can lead to a small but non-zero mass fermion. So we 



regard interactions ( |3.16|) and ( 3.17 ) as more plausible choices. 



The Higgs field is confined on the 3-brane. Its action is given by the brane action 

5Higgs = j d^x./^,[gZD^^^D,^ - Mo($^$ - ^^o')'} , (3-21) 

where = igWj[{x)^ + ig' B^{x) is the gauge covariant derivative. 

Now we have completed the model building in the 6D bulk. From the above, we know that 
the gauge field and Higgs field contents are the same with that in SM. The gauge-fixing terms 
and the ghost fields can be introduced as in SM. We omit them in this paper. 
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3.2 4D effective theories from the 6D bulk model 



We constructed the 6D bulk model in the last subsection. In this subsection, we plan to derive 



4D effective theories from the 6D bulk model in subsection 3.1 by KK decompositions. In 
that 6D bulk model, the gauge fields and the Higgs fields are confined on the brane; while the 
fermion fields propagate in the bulk. So we only need to reduce the fermion fields from 6D to 
4D. The process of reducing fermion fields and relative problems have been discussed in details 
in section ^. In this subsection, we derive 4D effective theories following discussions in section 
|2[ In subsection 3.2.1 , we give general results when we reduce the fermion actions from 6D to 



5D by KK decompositions. In subsection 3.2.2 , we discuss the special metric Eq. (|3.2| ). In this 



example, as we analyzed in subsection 2.3, we can obtain 3 families fermions in 5D by adjusting 



parameters in the metric. In subsection |3.2.3| , we further reduce the actions from 5D to 4D. By 
this step, we can obtain zero modes in 4D. These zero modes produce the 3 family fermions in 
SM by coupling with the Higgs fields on the brane. 



3.2.1 General discussions about KK decompositions from 6D to 5D 



In this subsection, we follow discussions in subsection 2.1. According to Eq. ( |2.4p , we expand 
the 6D fields U, D, U and V with 5D fields as 



2 I ^xvi^^^y^z) 



=5]G^(z)^/„(x^2/), (3.22) 



.2 



(x^ y,z) = Y, F2{z)Dn{x^, y), X^(x^ y, z) = ^ G^(z)D„(x^ y), (3.23) 



U 



V 



Xu 



2 I ^xui^^^y^^) 

Xu / 



^F^(z)^„(x^y),x^(x^y,z) = ^G^(z)^„(x^y), (3.24) 



^\ ) , x^(x^ y,z) = Y, F^{.z)&n[x^ y),xl{x^,y. ^) = E y)- (3-25) 

, / n n 



Note that we expand U and B with the same functions F,^{z) and G^{z). Because they are in 



the same doublet Q of SU(2)i, they have the same bulk mass parameters by Eq. (3.6). So they 
have the same expanding functions. These expanding modes are determined by equations 



^ - m^B^ F^{z) + (A* - M^)Gt{z) = 0, 
^ + m^B^ G^iz) - {X^ - M^)F^{z) = 0, 



(3.26) 
(3.27) 



where ^ can stand for Q, U or D, and we have employed the definitions 



F^iz) = B{z)-'/^F^{z), G*(.) = B{zr'/'G^J 



(3.28) 
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By these expanding, the 6D action (|3.6|) for Q is reduced to be 



ran / ^ \ ra n ' 



dz 



i ™ J- ^ ~l 



m n ' ^ m ^ n 



(3.29) 
, (3.30) 



where we have combined Un and Dn into a doublet Qn, because they have the same expanding 
functions. The action (3^) for U is reduced to be 



- / d'^xdy MH^A^i^rn'^n, 



Kin = Jdz [Ft*]^ + G^^Gl) , Mtn = j dz + G'i^G'i) 



(3.31) 



. (3.32) 



The action (^) for D is reduced to be 



Sv = j d\dy KZn 1^^^ [^m7'55^„ - 55#^7'^n + ^m7^5/.^n - 5^#™7^^n] 



i^^„ = j dz [F^*F^ + , M, 



(3.33) 
. (3.34) 



Now we consider the interaction sectors under the KK decompositions. For the gauge 
interaction sector ( p.ll[ ), ( p. 12 ) and ( p. 13 ), by the expanding in equations ( |3.22D -( 335| ), we 
obtain 



Suint 



9W^^{x)- + g'B^{x) 



d\dyA{yfKl^ |Q„7^ 
d^xdyA{yfK\i,^ {^n.l^{(j B ^{x)]^^] , 
d^xdyA{yfKZn {^^ml^[g' B ^{x)]%,] , 



Qn 



(3.35) 
(3.36) 
(3.37) 



where matrices K^, and have been defined by equations ( |3.3C1| ), ( ^.32 ) and ( 3.34 ) respec- 
tively. For the the Yukawa interaction sector ( |3.16D and ( 3.17| ), after the KK decompositions. 
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we obtain 



^OTukawa = j <fxdyA{yf \^^[y^n + i/?2?y^„]Q™^„$(x) ^<5(y - L') + H.C.j , (3.38) 

>^mn ~ j dzB{z) (^F^ + ' ^mn ~ j dzB{z) (^F^ F^ — , (3.39) 

^wYukawa = j d^xdyA{yf + ifhYi^^]Q^^nHx)^5{y - L') + H.C.j , (3.40) 

= / dzB{z) {F^^G'i + Gl*F^) ,Y^^ = J dzB{z) {f^*F^ - . (3.41) 

In the above, we obtain 5D effective actions by KK decompositions. We liave not consider 
the concrete form of the metric B{z). In next subsection, we choose the metric B{z) to be that 
in Eq. ( |3.2| ), and discuss the further simphfications of the above 5D effective actions. 



3.2.2 Deriving 5D effective theories for finite families 



In subsection [3.2.1 , we obtain the 5D effective actions for a general metric B(z). In this 
subsection, we plan to obtain 5D effective actions which include only finite KK modes. As 
we analyzed in section |2|, one can obtain finite KK modes by choose a special form of the 
metric B{z). In this paper, we choose the metric B{z) to be that in Eq. ( |3.2D . This metric 
has been analyzed in subsection 2.3 and in our previous paper |jT^. The results are that 
we can cut off the hypergeometrical series by requiring that it is finite when z ^ oo. This 
requirement determines the solutions uniquely up to the normalization constants. So we do 
not have freedom to imposing boundary conditions at the other boundary z = R. This implies 
that the normalization conditions Eq. ( 2.14| ) are not be satisfied, and we must change to the 



case (II) in subsection 2.1. In this case, we should redefine the fermion fields to obtain the 
conventional 5D effective actions like that in Eq. ( p. 20 ). As in subsection 2.1, we can make 
these redefinitions in two steps. 

Step (I): At this step, we analyze the kinetic terms of fermion actions. As the step (1) in 
subsection p.l| , we make the Cholesky decompositions for matrices K in the kinetic terms as 
follows. 



K 



u 



u 



K 



V 



V 



(3.42) 



One can make Cholesky decomposition for matrix K only when is a positive-definite hermi- 
tian matrix. In the numerical examples in subsection p. 3.1 and subsection 3.3.2, this condition 
is satisfied. Redefine the fermion fields as 



(3.43) 



In these new basis, the kinetic terms become the conventional ones similar to that of Eq. (|2.2 
while the mass terms are modified to be 



-1 



u-i 



MH 



v-i 



(3.44) 
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The fermions actions in the new basis are given by 



Sq 



d xdy 



(3.45) 



where Q can be Q, lA or P; while the corresponding ^ can stand for Q, ^ or ^ respectively. 

For the gauge interaction sector ( 3.35 ), ( |3.36|) and ( 3.37| ), by the definitions (l3.43| ), we 
obtain 



'S'Qint 
Suint 



d^xdyAiy)^ Q^^^ 



Qr 



d^xdyA{yf [^nl^[g'B^{x)]'Wn} , 
d^xdyA{yY [^nl^[g'B^{x)]^n} , 



(3.46) 
(3.47) 
(3.48) 



For the Yukawa interaction sector, after redefinitions (^.431), we obtain 



ukawa 



d'xdyA{yf { ^[y^^ + i(3vY^JQ^&nHx)^6{y - L') + H.C. 



L' 



where the index n is summed. Like matrices in the kinetic terms, the matrices in this sector 
become identities by the field redefinitions in Eq. ( |3.43D . 



(3.49) 
(3.50) 
(3.51) 
(3.52) 



^ L'-' 'mil ' "I- ' nirii rn-^ it - / 



SuY^ 



ukawa 



fu 



L' 



d\dyA{yf '-^\yln + iiiuyiLn\QA^{x)-rKy - L') + H.C. 



L'-' run • "i-iyi ' nirii -^^m - v— / 



The Yukawa interaction matrices are modified by the field redefinitions Eq. ( |3.43| ). 

In the above, we have completed the first step. This step is to make the kinetic terms of 
fermion actions to be the conventional ones. The mass matrices and the interaction sectors are 
modified accordingly. Especially, the gauge interaction sector becomes the flavor universal one, 
which is important to ensure the unitarity of the mixing matrix in the ZMA approach. 

Step (II): At the second step, we diagonalize the mass marix in the action ( p. 45 ). These 
matrices are hermitian, as they are defined in Eq. ( p. 44 ). They are diagonalized as 

(3.53) 
(3.54) 
(3.55) 



^ ^AqU^, AQ = diag(A^,A^,--- 

M" = U^^AuU^, Au = diag(A^, A^, • • • , A^), 



M 



V 



U^^AvU^, Ai, = diag(Af,A|', 



Redefining the fields in Eq. ( 3.45 ) as 



(3.56) 



These transformations are unitary. So the kinetic terms keep invariant; while the mass terms 
become the diagonal ones. By these transformations, the action ( p. 45 ) becomes the conventional 
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one 



(3.57) 



For the gauge interaction sector, by the unitary transformations (3.56), we obtain 



5'Qint 
5'Dint 



<fxdyA{yf [€nr[g'B^,{x)]€n] , 
d^xdyA{yf [§nl^[9B^{x)]§n] • 



Qr 



(3.58) 
(3.59) 
(3.60) 



Because the transformations in Eq. ( 3.56| ) are unitary. They keep universal for the flavors still. 
For the Yukawa interaction sector, after redefinitions Eq. ( ^.56 ), we obtain 



•^DYukawa 
yD 

•SwYukawa 
yU 



d^xdyA{yf { ^[5)^„ + iPvYZ]^^^nHx)^Ky - L') + H.C. \ , (3.61) 



fu 



D 



(3.62) 



L' 



d'xdyA{yf ^[y^^ + if3uY^JQ^^n^x)-6{y - L') + H.C. , (3.63) 



A 



L 



(3.64) 



The interaction matrices are modified by the unitary transformations ( 3.56 ). 

Now we complete the second step. This step makes the fermion mass terms to be diagonal 
ones. After this second step, we obtain the conventional 5D effective fermion action ( |3.57 ). 
The interaction sectors are modified by unitary transformations ( ^.56 ) accordingly. The gauge 
interaction sector are still universal for flavors after this step. 

We make some summaries about this subsection. By choosing the parameters in the metric, 
we can obtain finite KK modes. Because of this requirement, we must consider the normalization 



conditions case (II) in subsection 2.1. However, through twice redefinitions of fermion fields, we 
can also obtain the conventional 5D effective fermion actions. Having obtained the 4D effective 
actions for finite KK modes, we can derive 4D effective actions from these 5D actions in the 
next subsection. 



3.2.3 4D effective theories from 5D effective theories 

In this subsection, we further reduce the actions from 5D to 4D by KK decompositions. 
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We begin with the 5D effective actions obtained in last subsection. As that in subsection 



2.4, we expand the 5D fields with the 4D fields as follows 



Un{x,y) 


= E-4- 

j 




'f^^i}{y)4\{^) + &\y)^%{^) 


(3.65) 


Dnix,y) 


j 


-\y) 


\fil\y)d'ti^) + f^!i^{y)d%{x)_ 


(3.66) 


€nix,y) 


= E-4- 

j 


'\y) 


\fZi^\y)^f,L{-) + fZR^y^nU-) 


(3.67) 


^n{x,y) 




-\y) 


'f!f{y)d'ti-) + fnS^y)d%i-) 


(3.68) 



We give some interpretations about these expanding here. In the above expanding, the su- 
perscript j stands for different KK modes, while the subscript n can be interpreted as the 
family index. Note that they are not summed. We have expanded Un{x,y) and Dn{x,y) with 
the same functions fn^i\y) and f^'iiiy): because they have the same bulk mass parameters 
as shown in last subsection. As in subsection we require that these expanding functions 
satisfy equations 



dy 

d 



+ X^Ay:f{y)-^t^^)f:%\y) = ^^ 

^y->^^A^f:%\y) + ^l'^'f:fiy) = ^. 



(3.69) 
(3.70) 



where Q can be U or D as in last subsection; while ^ stands for Q, or ^ accordingly. 
Note that n can be regarded as the family index here and is not summed. For functions f^^'' (y) 
and fn^^{y): we designate the boundary conditions as in Eq. ( |2.62| ) 

+ A^a) f^^l^iy) = 0, /„«S^)(y) = 0, at y = L and L' . (3.71) 

These boundary conditions save left-handed zero modes while kill right-handed ones. These 
left-handed zero modes make the doublet of SU(2)i. For functions fZL{R)iy) ^^'^ fnL(R)^y^' 
we designate the boundary conditions as in Eq. ( p. 631 ) 



^ - ) CJ\y) = 0, fS\y) = 0, at y = L and L', 

d 
dy 



An>i)/S^(y) = 0, 4T(y) = 0, at y = L and L' . 



(3.72) 
(3.73) 



These boundary conditions save right-handed zero modes while kill left-handed ones. These 



right-handed zero modes make the singlets of SU(2)x,. As we discussed in subsection ^ 



these boundary conditions ensure that the expanding functions satisfy following normalization 
conditions 



£'dy[fTf:^^)=5-, £'dy[ 



dy (f:S'& 



(3.74) 
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where ^ stands for Q, or . Note that n is the family index and is not summed. 

By the above expanding ( 3.65| )-( |3.6§| ) and the normahzation conditions ( |3.74| ), the fermion 
action ( 3.57] ) becomes 



X \ - 
2 



■^(^■)7^a^^(^') - d^^^h^i^^^ - ^^i^^ , (3.75) 



where V' can be u or d. This action includes zero modes and massive modes. The modes 
and dn''^ are massless here. They obtain mass by coupling with Higgs field as in actions ( |3.61j ) 
and ( p3| ). 

For the gauge interaction sector, after the above expanding, we obtain 



5( 



Qint 



a 



Qn^ ^ 1 Qn,L 



Un,L 
dn,L 



Si 



Wint 



int 



(3.76) 
(3.77) 
(3.78) 



where n can be regarded as the family index and is summed. In the above, we have employed 
the normalization conditions ( |3.74 ). We have omitted three total minus signs in the above 
equations. We also isolate zero modes from massive modes obviously. The gauge interactions 
of zero modes are chiral because of the boundary conditions ( 3.71| ), ( 3.72 ) and ( 3.73 ); while the 
gauge interactions of massive modes are vector-like. We also see that the gauge interactions are 
universal for zero modes. 

For the Yukawa interaction sector, by the above expanding, we obtain 

a xayA—[y^^ + ipv^mn\2_^\Im,L Jn,R ^m,L^n,R + Jm,R Jn,L ^m,R"'n,Lj 



'CYukawa 



X <^{x)—6{y - L') +R.C. 

Ij 



(3.79) 



Si 



WYukawa 



^"■y^ j^i^mn + ^PU^mni 2_^\Jm,L Jn,R ^m,L^n,R + Jm,R Jn,L ^m,R^n,Lj 



X ^x)^6{y-L')+R.C., 



(3.80) 



where y^^^^ and y^(^) are defined as in equations ( |3.62| ) and ( p.64| ). Here m and n are the 
family indices. They are summed implicitly. These interaction terms include zero modes and 
massive modes. 

3.2.4 Mass matrices and mixing matrix 

In the last subsection, we have derived the 4D effective actions from the 5D ones in subsection 
3.2.2| . In this subsection, we derive the mass matrix for 4D zero modes and their mixing matrix. 
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Before doing that, we convert the gauge field action and the Higgs field action to the canonical 
forms. 

For the gauge field, use the metric ( |3.4| ), the action ( 3.14| ) becomes to be 



Sgauge — I d X < Tj^ Tj ^ 



(3.81) 



where ry'^" = diag(— 1, 1, 1, 1) is the 4D Lorentz metric. Note that we do not need to redefine 
the gauge fields. So the gauge interaction actions ( |3.76| ), ( |3.77| ) and ( 3.78| ) keep invariant and 
still apply in this subsection. 

In order to convert the Higgs field action to the canonical form, we redefine the Higgs field 

as 

^(x) = ci>(,), ~s = s'^^^. (3.82) 

By this redefinition and the metric (p.4|), the Higgs action ( 3.21| ) changes to 

^Higgs = j d^x [r]'"'D^^^D,ip - fioi^^if - vY] , (3.83) 

where v = s^vq. For sjj <C 1, it supplies a beautiful geometrical solution for gauge hierarchy 
problem suggested by Randall and Sundrum in Because the gauge fields do not need to 
be redefined, the gauge covariance derivative keeps with the same form as that in ( |3.21| ). 

After the redefinition ( p. 82] ) for Higgs filed, the Yukawa interaction terms ( |3.79| ) and ( p.80| ) 
change to be 

•^VYnhnwa " / " ^ JA " ^PT^^mnl / , \ Jm.L Jn,R ^m,L"'n,R + Jm,R Jn,L ^m,R"-n,Lj 

X ^{x) (^^) 5(y-L') + H.C., (3.84) 

'!5WYukawa " / « XdyA—[y^^ + lljuymn\2_^\Jm,L Jn,R Qm,L^n,R + Jm,R Jn,L Qm,R^n,Lj 

x^{x)(^^^ 6{y - L') + R.C.. (3.85) 

Following the ZMA approach, we isolate the zero mode terms from above expressions as 
follows, 

+ H.C., (3.86) 

Wawa = / d^-dyA^[yl!^^ + ^PuY^J{fSy:PQt^^^^^^^ (t)''^^"^') 

+ H.C.. (3.87) 

As in SM, after that the Higgs filed develops a vacuum expectation value, the electroweak 
symmetry breaks. These Yukawa interactions produce mass terms for fermions. From the 
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above, we see that the mass terms are related to the fermion zero mode profiles. We have 



worked out these profiles in subsection 2.4 and give their approximation behavior there. By 



these profiles, we obtain the mass matrices for quarks as follows 



M' 



d 



mn 



-^j^lJ^mn ~^ ''PVJ:mnlJm,L Jn,R ' 1"^' 



= + (3.89) 

where the indices m and n are not summed. The term Ymnfmfn stands for the product of three 
quantities Ymn, fm and /„. v is the Higgs vacuum expectation value as in Eq. ( 3.83| ). fn^i\ 

fnR^ fnR^ valucs of the canonical zero mode profiles as we defined in equations 

( |2.67 ) and ( 2.6^ ) in subsection 2.4. They are given by 



fS = ^^^^/--''^ c^ = A^L, (3.90) 



?^(0) _ / -2cf-l cf-l -\^T C^QU 

Jn,R ~ \l i_ g-2ef-l ^ ' ^n-K^: (,3.9ij 



Jn,R - Y i_g-2e-l' (3.92) 



where A^, and A^ are determined by equations (p.53|) , (|3.54|) and ( ^.55 ). We can rewrite 
above equations with the matrix form as 

= diag(/^^r , . . . , fS), Pi = diag(/?f , • • • , /3), (3.94) 
= ^PLiy^ + iPuY^]PR, (3.95) 
= diag(/?f , . . . , fS), Pi = diag(/S^ • • • , /3). (3.96) 

These mass matrices are general complex matrices, and are not hermitian matrices. We can 
make the single-value decompositions for them to derive the mass eigenstates as follows 

M'^ = vl^MdVdR, Md = diag(md,,md2, • • • (3.97) 
M" = vl^MuVuR, Mu = diag(m„,,m„2, • • • ,m„„), (3.98) 

where md^(^Ui) > for i = 1, 2, • • • , n according to the definition of single-value decomposition. 

As we discussed above, the gauge interaction terms keep the same form with that of equa- 
tions ( 3.76| ), ( |3.77|) and ( ^.78 ). The zero modes interact with gauge fields just like that in SM. 



So we can define the mixing matrix for quarks like that in SM as 

VcKM = VulVJ^. (3.99) 
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3.3 Numerical results 



In subsection 3.1, we construct our model in 6D bulk. In subsection we derive 4D effective 
actions from the 6D ones by two step KK decompositions. In this subsection, we give numerical 
examples to show that the results of our model can be very close to the experimental data. 



Prom the model in subsection |3.l| , we know that there are many parameters in this model. 
These parameters are not determined by the model. We need to input these parameters by 
hand to obtain the numerical results. These parameters can be classified to two groups: the 
parameters in the metric and the parameters closely related to the fermion mass matrices. 
Before giving the numerical results, we give some discussions about the permissible extent of 
the parameters. 

We discuss the parameters in the metric at first. For the extent of the extra dimension y, 
we let 

e = ^ = 10~^^ y = 10^9 GeV, ^ = TeV, (3.100) 



which is necessary to interpret the gauge hierarchy as suggested in |13]. From Eq. (3.82), we 



know that vq relates to v by the factor sjj. The choice of Eq. ( 3.100| ) implies that s ~ 1 in 



order to interpret the gauge hierarchy. The value of e coincides with that in the paper |]18| . We 
designate the boundary value of the dimension z by the equation 



40. (3.101) 



b 

We also designate the value of the parameter to in the metric B{z) by the equation 

ujL = 0.15. (3.102) 



Lo can be regarded as the intrinsic sale of the the dimension z. Eq. ( 3.102| ) implies that u; is 



about 10 percents of the Planck scale. We designate o in the metric B(z) by the equation 

^ = 6. (3.103) 

The other parameters like s and b are not necessary designated in the numerical examples as 
they always emerge in combinations with other parameters. 

In the above, we have designated the necessary parameters in the metric. We notice that 
the values of these parameters are not determined by our model. We expect that they can be 
determined by some underlying theories, which are not discussed by the present paper. We 
choose them to the above values by hand in this paper, because we find that these values can 
make the results of our model to be very close to the experimental data. 

3.3.1 Numerical results in quark sector 

In this subsection, we give numerical results in the quark sector. The parameters in the metric 
have been given above. To obtain the numerical results, we need to further designate the 
parameters related to quark mass matrices. At first, we need to designate the parameters m 



and Ai in actions ( |3.6D , ( |3.7D and (p^). Prom the analysis in subsection 2^, we know that the 
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Table. 3.3.1 Parameters in quark sector 



Field 

q 0.152 4.25 

U -0.173 -4.0 

V -0.093 -4.2 

/w=15.5 TeV (5u = 0.55 - 0.43z 

fv=l'2.7 TeV Pt) = -1-2 - 0.5i 



number of family is determined by the pair (f , ^s)- So the parameter is closely related to 
the number of family. In subsection we have given a parameter set which permits the very 
3 families. As we analyzed in subsection |2.5| , these conclusions also apply to the new setup in 
subsection 'LB. This parameter set is given by equations ( ^.55 ), ( 2.56| ) and ( p. 57 ). We have 
designated the value of | by Eq. ( 3.103|) , so the possible extent for is given by 



m 

—s e 



1 



35 + 5' + 



and —s + 

10 



(3.104) 



We exclude the point ^rs 



10' 



because when 



j_ 

10' 



we obtain a massless solution, which 



coincides with the zero mode solution as we discussed in subsection |2.3|. The parameters m in 



actions (3.6), ( |3.7| ) and (|3.8|) should take values in the intervals Eq. (3.104) to ensure that there 
are the very 3 families in our model. While M. is irrelevant to the family number. Instead it is 
closely relevant to the values of quark masses, as it can be seen from the following numerical 
example. For ra in the intervals Eq. ( |3.104| ), the explicit expressions for 3 family KK modes can 



be determined as we discussed in subsection 2.5. According to those discussions, the solutions 
are very similar to that given in appendix The differences are that we should replace A 
with A A in the expressions in appendix p. By the normalization conditions Eq. ( B.25 ), these 
solutions can be determined completely. The values of m and ^A for different fields are different 
generally. We adjust them by hand to fit the experimental data. We give the values for them in 



Table. 3.3.1. Note that the negative real numbers m also emerge in Table. 3.3.1. The solutions 
for this case are also discussed in appendix ^ 

The Yukawa couplings / and /3 in equations ( ^.16 ) and ( 3.17 ) also need to be input by 



hand. We give their values in Table. 3.3.1. In Table. 3.3.1, we have defined that 



AL e'^^ + a 



AL e"^^ + a 



(3.105) 



Having designated these parameters, we can obtain the numerical expressions of kinds of 
quantities in our model, like the matrices K and mass matrices M in the fermion actions ( |3.29| ), 
( 3.31 ) and (^3.33 ), the eigenvalues A„ in equation (|3.57| ) after two step redefinitions of fermion 
fields and so on. In this paper, we omit these intermediate numerical expressions. We only give 
the final mass matrices in equations ( 0.93 ) and ( p. 95 ). However, some analytical expressions 



for these intermediate quantities can be found in subsection 4.1 
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By the parameter values given above, we obtain the numerical expressions for mass matrices 
dPD and dpI) as follows. 



/ (6.44 + 0.552i) • 10-"^ (0.922 + lA9i) ■ lO'^ -0.123 - 0.158i 
(-1.53 - 3.25i) • 10~3 0.573 -5.44 - 4.98i 

(1.29 + 1.64f) • 10-2 2.23 + 2. 36i 135.58 - 13.4i 

/ (0.72 + 2.78i) • 10-3 (1.1 -0.58i) • 10-2 (3.9 - 7.85i) • lO'^ 

(-6.57 + 1.8i) • 10-3 4.39-10-2 (-1.36 - 0.49i) • 10-^ 
y 0.155 -0.46z -0.984 -0.93i -1.30 + 1.19? 



GeV, (3.106) 



GeV. (3.107) 



Making the single- value decompositions as in equations ( 3.97 ) and ( ^3.981) , we obtain the quark 
masses as follows 



m„ = 0.64 MeV, = 584.87 MeV, mt = 136.48 GeV, 
md = 2 MeV, = 36.36 MeV, rrib = 2.278 GeV. 



(3.108) 
(3.109) 



They are consistent with the MS quark masses evaluated at 1.5 TeV in the paper |18]. By 



Eq. (3.99), the mixing matrix and its absolute value are given by 

VcKM 



VcKM 



/ 0.974 + 0.00614i 0.224545 - 0.0293i (4.297 - 0.58i) • 10-^ 

0.1488 - 0.1699i -0.535961 + 0.811i 0.0553 - 0.0232i 
y -0.0011 -0.0172i 0.01957 + 0.0542i -0.249 + 0.9666i 

0.974014 0.226448 0.00433647" 
0.225832 0.972317 0.0599961 
0.0172432 0.0576282 0.998189 



(3.110) 



(3.111) 



We also obtain the Jarlskog invariant as 

J = -lm{VudVcbV:f,V*a) = 3.20416 x IQ-^ 



(3.112) 



They are very close to the experimental data compiled in [19| 



3.3.2 Numerical results in lepton sector 

In this subsection, we discuss the lepton sector. We suppose that neutrinos are Dirac ones. In 
this case, the lepton sector is very similar to the quark sector. For neutrinos in other scenarios, 



see ||0|. 

We introduce the fermion field contents in the 6D bulk as 



L 




M, £. 



(3.113) 



They transform under the gauge group SU(3)c x SU(2)j^ x U(l)y as L = (l,2)_i, M = 
(l,l)o, £ = (1,1)_2. Note that we also introduce only 1 family lepton in the 6D bulk. The 
actions of these fields are the same with that in equations (|3.6D , (|3.7| ) and (^^) . The model in 



subsection 3.1 applies similarly here, other than there is no color interaction for leptons. The 
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Table. 3.3.2 Parameters in lepton sector 



Field 



M . 



L 
M 
S 

/V=23.89 TeV 
fe=2l.22, TeV 



0.16 -9.0 
-0.0944 -9.75 
-0.147 -5.6 
f3j^ = 1.5 + 0.45i 
(3e = 0.9752 



process of deriving 4D effective actions from the 6D ones in subsection also applies here. 
While the mixing matrix for leptons is defined by 



^PMNS = ^lK/L- 

Now we discuss the parameters in the lepton sector. The parameters in the metric given by 
equations ( 3.100| ), ( |3.101 ), ( 3.102 ) and ( 3.103 ) still apply in the lepton sector. The parameters 
m of leptons should also in the intervals Eq. (3.104) to ensure that we can obtain the very 3 
families in the lepton sector. Other parameters should also be input by hand like that in the 
quark sector. We adjust them by hand to fit the experimental data. We give their values in 
Table. |3.3.2 . Having designating these parameters in lepton sector, we can obtain the numerical 
expressions for kinds of quantities as in the quark sector. In this subsection, we also only give 
the mass matrices for leptons. They are given by 

/ 3.204 - 8.405i 21.18 - 13.78i 6.37 + 26.6i \ 
0.6978 + 0.632i -6.286 3.53 - 2.072i • 10^^ eV, (3.115) 
-6.48 + 7.72i -12.29 - U.lli 5.736 - 26.22ii j 

( 0.287455 + 3.48i -10.01 + 90.3488i 1501.6 + 154.77i \ 

-0.4266 - 6. 935ii 129.213 58.594 - 952.495i MeV. (3.116) 
10.40 - 1.072ii 9.527 + 85. 96i 37.56 - 454.82i / 



(3.114) 



v 



Making single-value decompositions for these matrices as in equations ( p. 97 ) and (3.^ 
obtain the lepton masses as follows 



we 



me 



0.511 MeV, 
0.0019 eV, 



ra,. 



7712 



= 105.229 MeV, 
0.013 eV, mg = 



rrir = 1849.15 MeV, 
: 0.05 eV. 



(3.117) 
(3.118) 



The masses of electron and muon are close to their experimental value; while the mass of r 
is moderately large than its experimental value compiled in [19|. The neutrino masses are of 
normal hierarchy type. They are close to the experimental values in P]. We can also obtain 



the mixing matrix defined in Eq. ( 3.114 ) and its absolute value as 

/ 0.6799 - 0.4592i 0.5163 0.2442i 0.02504 0.003909A 

-0.4619 - 0.02652i 0.2922 + 0.5664i 0.44075 + 0.43071i , (3.119) 
0.2051 - 0.2659i -0.5136 - 0.06183i 0.7715 - 0.156H / 

/ 0.820465 0.571134 0.02534 \ 
0.462666 0.637319 0.61625 . (3.120) 
0.335822 0.517329 0.78714 / 



^PMNS 



VpMNS 



V 
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They are in the 3a extent of the experimental values in |^. 

We have supposed neutrinos to be Dirac ones. So as in the quark sector, we can calculate 
the Jarlskog invariant for the mixing matrix as 

J = -Im(KiV;.3K3^;i) = 19.193 X 10-^ (3.121) 

It is lager than that in the quark sector. We note that its size is not the inevitable result of 
our model. When we adjust the parameters, we find that the size of J is very sensitive to the 
mass m of the field M. By adjusting ^^s, we can still keep the mixing matrix and the neutrino 
masses to be close to their experimental values, but J can varies remarkably. 



3.3.3 Brief comments on the numerical results 

We give the numerical results in the last two subsections. We see that there are still many 
parameters in our model. The parameters are closely relevant to the absolute size of the 
fermion masses. The parameter of the field M is more lager than that of quarks and charged 
leptons, because the neutrino masses are remarkably smaller than the masses of quarks and 
charged leptons. While the parameters are closely relevant to the mass hierarchy structure 
of quarks and charged leptons. In general, lager absolute value of produces larger mass 
hierarchy. The parameters P are closely relative to the CP violation measure J. /3 = implies 
that J vanishes. In the above numerical examples, we see that /3 7^ in the lepton sector, so 
CP violation also emerges in the lepton sector. We have not found a group of parameters with 
/? = which can fit the experimental data as the group of parameters in the last subsection. 

We also see that the values of —s are about 100 times larger than the values of —s, while 
this little hierarchy is not explained in this paper. In addition, because there are too many 
parameters in our model, we only adjust them by hand to obtain the numerical results. These 
numerical results are very close to the experimental values. We have not done further adjust- 
ments to make them in the la extent permitted by experiments. In the last two subsections, we 
only give rough numerical examples to show that our model can close to the experimental data 
in high precision. For the quark masses, we take the renormalization effects into consideration, 
and adjust the parameters to fit the running masses at 1.5 TeV compiled in [p!^]. While for 
the leptons masses and the mixing matrices, we only adjust parameters by hand to fit the face 
values compiled in ^ |l^, and the renormalization effects are omitted. So the above numerical 
examples are only rough treatments. The renormalization effects of these quantities should be 
considered for more detailed comparison with the experimental data. 



4. Some analytical treatments about the model and more relevant discussions 

In section ^, we introduce our model and give the numerical results. At first sight, this model 
seems very complicated. However, when we give some analytical treatments about this model, 
we will see that some quantities in this model have very concise expressions. We can make some 
qualitative conclusions from these concise expressions. These analytical treatments can help us 
to understand how our model works more clearly. In subsection O, we make more discussions 
about several relevant problems. 
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4.1 Some analytical treatments about the model 



For analytical discussions, we take fields Q and T) for example. The analytical treatments for lA 
are very similar to that of Q and P. As we discussed in appendix^ the solutions for KK modes 
are different for m to be positive or negative. While the eigenvalues A determined by equations 
(I2I9I) , (|2^) , (|2^ and (|2^ ) can be real or pure imaginary according to the parameter pair 
(|, ^s). The analytical expressions can be classified into four classes: (1) m > 0, A is real; 
(2) m > 0, A is pure imaginary; (3) m < 0, A is real; (4) m < 0, A is pure imaginary. For the 



new setup in subsection |2.5| , we should replace A with A A as we analyzed in subsection p.5| . In 
the numerical examples for quark sector in subsection p.3.1| . The solutions for Q belong to the 
class (2); while the solutions for V belong to the class (3). The solutions for U belong to the 
class (4). In the following discussions, we always suppose that the solutions for Q belong to 
the class (2) and the solutions for V belong to the class (3). Other situations can be discussed 
similarly. 

For the analytical treatments, we must discuss the solutions for KK modes at first. The 
solutions can be determined according to our discussions in subsection 2.5. They have the 
similar forms and characters to that given in appendix ^. By these solutions, we obtain the 
expressions for the matrices K and M in the fermion actions ( 3.29] ) and ( |3.33| ) as follows 
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= MvK^ 
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AA^ 









\av 





-AA^ 



(4.1) 



(4.2) 



where we have used the normalization conditions in Eq. ( p.25| ). We give the expressions for 
oq, a-Di bq and 6d in equations ( p.l| ) and ( p.2| ) in appendix |y. They are all real numbers 
by definitions. Here and in the following, we arrange the column and the row indices for the 
matrices as 0, 1, —1, which are the indices for KK modes displayed in appendix M. AA are the 



eigenvalues in equations ( 2.79 ) and ( 2.80 ). They can be calculated according to our discussions 



in subsection 2^. We have employed Eq. ( 2.78 ) to rewrite the expressions for mass matrices 
M. 

Now we discuss the matrices in interaction sectors. For the gauge interaction sector, from 
equations ( |3.35| ) and ( ^.37 ), we know that the matrices in this sector are the same with that 
in the above. The matrices in the Yukawa interaction sector are important for our discussions. 
We obtain expressions for these matrices in Eq. ( 3.39| ) as follows 
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(4.3) 



The elements in these matrices are defined by equations 
elements can be complex numbers generically. 



and ( C.4 ) in appendix^. These 
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Prom the analytical expressions above, we see that these matrices all have very concise 
structures. Following the procedures in subsection 3.2.2, we need two steps to obtain conven- 
tional 5D effective actions. 

Step (I): Making Cholesky decompositions for the matrices K, for the field Q, we obtain 



diag(A^,A^,A3 







-i 0^ 
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ly 






(4.4) 
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+ 6 



While for the field P, similarly we obtain 
= diag(Af ,A 



A^) 



cos t/x) 

1 



V2 



sm Ux> 



72 """'^^ 75, 



V 



= diag 
, tan 9x> 



(4.5) 



av + 



^^2 1 



8a 



V 



+ 6.2 



The expressions for these eigenvalues in the above are given by equations ( p.5D and ( p.6| ) in 
appendix y. By the redefinitions of fermion fields in Eq. ( |3.43 ), we know that the matrices 
in the kinetic terms of fermion actions and gauge interaction terms become to be the identity 
matrices. While the matrices in the mass terms become to be 



Mn + 







(4.6) 



Mv + 



(4.7) 



-ibi 
y -ioQ iftdj 

/ a©' 
6© 

The elements of these matrices are defined by equations ( p.7|) and (|C.8| ) in appendix |^. By 
definitions, aq, bq, ax) and bx> are all real numbers. We see that the above expressions are still 
of very concise structure after this first field redefinitions. We do not give explicit expressions 
for the Yukawa interaction sector here. We will give them in step (II). 

Step (II): As in equations (|3.53| ) and (|3.55| ), we diagonalizing the matrices in equations 
(14. 6|) and (14.71). For the mass matrix of field O, we obtain 
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(4.8) 



^ sin i^fT 



'V2 
1 

'V2 



tan I 



1 

V2. 
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While for the mass matrix of field T>, we obtain 



= U^^AvU^, Ax, = diag(Af , A^, • 
/ cos-dv -jK sin tDti — j^s'mi^Ti' 



sin 



V2- 
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71 
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V2' 
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V2 



COS ■& 
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75^ 



AA^, Xf 



bv 



Mv + ^X^, aP 



1,^ 



+ 62 



(4.9) 



These matrices are still of concise structure. If = 0, we see that the eigenvalues A have the 



spectrum 0, ibAA. This spectrum is similar to Eq. ( 2.45 ), which is the spectrum of A before the 
redefinitions of fermion fields. While for M ^ 0, the spectrum is similar to Eq. (2.82). These 
results are expected in subsection |2.5| . 

By the redefinitions of fermion fields as in Eq. ( p. 56 ), we know that the fermion actions 
become to be the conventional ones as in Eq. ( p. 57] ); while the gauge interaction terms keep 
to be the flavor universal ones as in equations (|3.58| ) and ( |3.6Cl| ). For the Yukawa interaction 
sector, after the redefinitions, we obtain final results of these matrices as follows 
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(4.10) 
(4.11) 
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where 



and 



are the inverse matrices of 



and V A-^ respectively. We define the 



elements of matrices H'^ and in equations (Icl ), (IC.IOD , (|C.11| ) and (|C.12D . From the 
above, we see that the transformation matrices H'^ and have the concise structure similar 
to that of matrices and in Eq. (|4.3|). Working out the product of these matrices in 



Eq. ( 4.10 ), we obtain 
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(4.13) 



The elements in these matrices can be expressed with the quantities in the matrices H'^, , 
y^ and Y^ by the matrix multiplication. We omit their explicit expressions for simplicity. 
Because of the special structure of the transformation matrices H'^ and , the structure of 
y^ and Y^ keep invariant under these transformations and only their elements are modified 
to be different values. 
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Making single- value decompositions for matrices and , for , we obtain 
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sin dy 



sin 6y 



diag(E^, S^, S^), 
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For Y^, similarly we obtain 
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Y 



Vyl 




s^=diag(Er, sr, 

sin 6y \ 



Sr), (4.15) 



7S"°'^^ 7S 



Some quantities in these expressions are defined in equations ( |C.13 ), (|C.14 ) and ( C.15 ) in 
appendix The unitary matrices Vyn and Vyr have the similar structure to that ofVyi and 
Vyl, we do not give them explicitly here. 

In the above, we give some analytical treatment about our model. We only display the 
results for down quark sector, but the results for up quark sector are similar to that in the 



down quark sector. As we discussed in subsection 3.2.2 , the lepton sector is similar to the quark 
sector. So the above discussions also apply to the lepton sector. From the above discussions, we 
see that the matrices in the 5D effective actions all have very concise structure. These concise 
structures are induced by the special characters of KK modes as we analyzed in subsection 2.3 
and subsection |2.5| . Especially, the unitary transformation matrices VyL and Vyl in equations 
(4.14) and ( 4.15| ) are very close to the structure of the experimental PMNS mixing matrix in 
the lepton sector. So the Yukawa coupling matrices y^ and Y^ are appropriate to construct 
models for lepton mixing matrices. However, the addition of these two matrices y^ + i[3x>Y^ 
in Eq. ( ^.61 ) becomes to be a general matrix. It does not have the concise structure like that 
of matrices and Y^ . Its eigenvectors are complicated and we do not give them here. 

The matrices y^ and Y^ are the 5D Yukawa couplings. The physical 4D mass matrix is 
given by Eq. (|^. From Eq. (|3^ ), we see that these matrices are modified further by the 
4D zero modes profiles. The concise structure of y^ and Y^ are lost and are distorted further 
by these zero mode profiles to be a general matrix. 

By the above analytical treatments, we can make some qualitative discussions about how 
our model works. From Eq. ( f4.3D , we know that the 5D effective Yukawa couplings are deter- 
mined by the profiles of KK modes. Due to the special characters of KK modes as we analyzed 
in subsection |2.3| and subsection 2.5, they have very concise structures. After two step field 



redefinitions, the induced Yukawa couplings Eq. ( 4.13| ) are still of concise structure. These 
concise structures are distorted by the summation in Eq. ( p. 61 ). When we reduce further the 
actions from 5D to 4D, the 4D zero mode profiles distort the 5d Yukawa couplings further. The 
exponential behaviors of 4D zero mode profiles also induce the hierarchy mass structure in 4D. 
We see that the structures of mixing matrices VyL and Vyl are universal for quark fields. In our 
model, because the lepton sector is similar to the quark sector, these structures also apply to 
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the lepton sector. These concise matrix structures are distorted by the two sources we discussed 
above: the summation + i(3x>Y^] and the 4D zero mode profiles in Eq. ( 3.93| ). These two 



sources distort these concise matrices to some general matrices. Their analytical expressions 
are complicated and are not appropriate to make qualitative discussions. The above discussions 
give some sketchy interpretations about how our model works. 

4.2 More relevant discussions 



In subsection 2.5, we suggest a new action Eq. ( 2.75| ) for our model building. This new action 



breaks the 6D local Lorentz invariance obviously. In this subsection, we discuss the possible 
origin of this action. 

The new term ^T'^'e^F'^^' suggests that we may introduce the gauge interaction term 



^7v*r'e;;r''^', (4.16) 

where Am is an Abelian gauge field. After that Am acquires the background value Aq = v, Ai = 
0, i = 1, 2, • • • ,5, this interaction term supplies the term v^T^ e^^T"^'^ . However, we argue that 



this interaction term is not a proper choice. Because T'' emerges in Eq. ( 4.16| ), we need to define 
the local gauge transformation for ^ as 



^^e'"' ^. (4.17) 

However, under such a gauge transformation, the term transforms as — > — ^f^. The 
mass terms is prohibited by this gauge transformation. So the interaction term Eq. ( [4.16| ) 
is not appropriate to produce the new action Eq. ( 2.751) . 



Instead of the above, we suggest another way to produce the action Eq. (|2.75|) . We introduce 
a 5-form gauge field strength as follows 

= -^FmnhklF^''''''^ (4.18) 
2 ■ 5! 

FmNHKL = d[M^NHKL]- 

This 5-form gauge field strength is similar to the Maxwell electromagnetic field strength. We 
suppose that this 5-form field interacts with the fermions by the nontrivial interaction term 

■^-t V^^^^MiVHi^Lr^^r^r^r^r^V', (4.i9) 

7^2 1 7T1 jpMNHKL 

r = —-^—^.tMNHKLr , 

where a is a constant real number. We have defined that = e^F*^, in which F*^ are given by 
Eq. ( |2.3| ) and the index a is summed. 

In order to make the 5-form field to produce the appropriate background value, we construct 
the following interaction system 

S = J d'^xdydz^/^{2M^R} + j d^xdydz,/^ !^^g^^VK^^L(t> + V{(t>)^ 

+ f d^xdydzV^l-^FMNHKLF^'''''''''] . (4.20) 
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in which (p is a scalar field and V{(j)) is its potential term. We solve this system supposing the 
metric Ansatz Eq. ( |2.lD . Suppose that 



f{z)e^-^P\ f{z)=vB~\z), 



-g 



(4.21) 



l2,i',a,P,-f = 0, 1,2,3,5, e 



01235 



1, 



where f is a constant. Then the equation of motion of this 5-form field strength and the Bianchi 
identities for it are both satisfied. We also suppose that the scalar field (p only depends on the 
coordinate z. By the metric Ansatz Eq. ( |2.lD , we obtain the following equations 
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(4.22) 
(4.23) 



+ (4-24) 



8M4 



(4.25) 



in which Ay = fi^ = (/)^ = ^. Obviously, A(y) should be of the form 



1 



ky + c 



(4.26) 



in which k, c are constants. This system is similar to that in our previous paper |1^. The 
conclusions for that system apply here. For any B{z), there exists an appropriate V{(j)), which 
makes Eqs. (|^)-(|^) to be satisfied. So we can choose an appropriate V{(l)) to make the 
metric Eq. ( |2.46| ) to be our background solutions. This system is also complicated, and we can 
not find a superpotential to express the general solutions. 

By the above discussions, if we suppose a = | in Eq. ( [4.19 ), then we obtain 



(4.27) 



This term is exactly that we expect in Eq. ( 2.75| ). So the above system can give the new action 
Eq. ( ^.751) and give the background solutions Eq. (|3.2D at the same time. However, note that the 
term Eq. ( [4.19| ) has a very nontrivial form for a = |, which implies that it will be intractable 
when we treat it as a quantum theory. We have not found a more simple term to replace it yet. 



4.3 Gauge fields in the bulk and a spontaneously broken framework for CP violation 

In our model building in subsection we have supposed that the gauge fields are confined on 
the 3-brane sited at {z = R, y = L'). By this assumption, the gauge fields only propagate in the 
physical 4D spacetime, and we do not need to make the conventional KK decompositions for 
them. In this subsection, we consider the possibility that the gauge fields propagate in higher 
dimensions. 
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When we consider gauge fields in higher dimensions, a natural choice is that gauge fields 
propagate in the 6D bulk. However, gauge fields propagating in the 6D bulk induces intractable 
problems. Because the 5th space dimension z is intrinsically semi- infinite in the metric Ansatz 
( |3.1| ), the bounded KK modes of gauge fields must be non-constant. Such non-constant profiles 
of gauge fields would break the unitarity of mixing matrix in the ZMA approach. This situation 
differs from the conventional flavor models in RS spacetime. In RS flavor models like P, 18 1, 



the 5th dimension is a finite interval and the zero modes of gauge fields have constant profiles. 
These constant zero modes profiles keep the unitarity of mixing matrix in the ZMA approach. 



So in our model building in subsection |3.1| , we consider the situation that the gauge fields only 
propagate in the 4D spacetime. This choice makes the numerical results of our model to be 
close to the experimental data. 

As we just discussed above, gauge flelds in the 6D bulk induce intractable problems. How- 
ever, the gauge flelds can propagate in the 5D spacetime. We can consider the situation that 
the gauge fields are confined on the 4-brane sited at z = i?. Because the zero mode proflles of 
gauge fields can be constant, the unitarity of mixing matrix can be kept in the ZMA approach. 
By some modifications, the model in subsection can apply similarly in this situation. 

Moreover, the gauge flelds propagating in 5D can bring interesting influence on the me- 



chanics for CP violation. In our model in subsection 3.1, CP violations originate from the terms 
after the coefficients /3, which are put in by hand. If we consider that gauge flelds propagate 
in 5D, we may have a new mechanics for CP violation. According to the mechanics suggested 



in 1 21], the 5th component of gauge field can develop a vacuum expectation value through the 



gauge invariant line integral 

{Ay) = j dyAy. (4.28) 

This vacuum expectation value can break the non-Abelian gauge group and also supply an 
origin for CP violation. In order to make a realistic model, we may embed the electroweak 
unification group SU{2) x U{1) into a larger unification group SO{5) x U{1). The vacuum 
expectation value in Eq. (|4.28D has also been discussed in gauge-Higgs unification framework 



22f| . So instead of putting in CP violation by hand as we did in subsection |3.l| , the gauge fields 
in 5D can supply a spontaneously broken framework for CP violation according to the above 
mechanics. 



5. Conclusions 

In warped extra dimensional RS model, the fermion mass hierarchies can be produced by the 5D 



bulk mass parameters of the same order. In our previous paper [13], we suggest that these 5D 
mass parameters can be interpreted in a two-layer warped 6D model, and such an approach also 
supply a solution for family problem. In this paper, we combine these suggestions and construct 
a specific model to address the fermion mass hierarchy problem and the family problems at the 
same time. We give numerical examples in subsection |3.3| to show that the numerical results 
of this model can be very close to the experimental data in both the quark sector and the 
lepton sector. However, because there still exist many parameters in our model, we only make 
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rough numerical treatments about the model, and do not further adjust parameters to fit the 
experimental data in higher precision. 



We further make some analytical treatments for our model in subsection 4.1. These an- 
alytical treatments show that some very concise structures exist in this model. They imply 
some common features shared by quarks and leptons. However, the breaking of these concise 
structures makes the model to be a complicated one, and we do not make more analytical dis- 
cussions. Some approximate treatments may be helpful to illuminate this model more clearly. 



In addition, a natural question is that whether we can interpret the parameters in Table. 3.3.1 



and Table. 3.3.2. It seems that it is appropriate to interpret the origin of those parameters in 
a grand unification framework like in [^3| . 
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A. Analysis about the equations determining the eigenvalues 



In this appendix, we analyze how equations ( |2.49D , ( |2.50| ), ( |2.51| ) and (2.52) restrict the number 
of eigenvalues to be finite. 

Solving equations (|2l9| ), (^ISOD , (|2]5|) and (^ ) for A, we obtain 
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for d^OD, 
for (pT|), 



for (12:521) . 



(A.l) 
(A.2) 
(A.3) 
(A.4) 



From these solutions, because we suppose that ^s, ^ are real and n is a non-negative integer, 
we can infer that (^)^ must be real, in other words, ^ must be real or pure imaginary. For 



simplicity, define t 



•A^2 



, then t is a real number. 



In the following analysis, we suppose that m > 0. For the case m < 0, we discuss it in 
appendix ^. Now we discuss these solutions in two cases: 

Case (1): In this case, let f > 1, so p = ^s(l - f ) < 0. Rewrite equations ( ^^491 ), (pISOl ), 
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(12^511) and ([2l5^ ) as 
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Because /O < and n > 0, Eq. ( A.8 ) has no solutions for any t obviously; while equations (| 

v2 



[A.q ) and ( |A.7| ) might have solutions only when ^ < (^s) • Define 
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f{t) and ^(i) are both increasing functions about t when t < (^s)^. By this feature, we obtain 

-1, (A.ll) 



< fit) < -s 

UJ 



m / /a\ 2 



By equations (A.ll) and (A. 12), we can determine the extent of n as 
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(A.12) 

(A.13) 
(A.14) 
(A.15) 



When f > 1, 

' Lu 



i + (f-i) 



< 0. So Eq. ( [A.15 ) is impossible, and the corre- 
sponding equation ( |A.7| ) has no solutions. 

In summaries, in the case, the extent of n is bounded by equations ( |A.13 ) and ( |A.14 ). Note 
that it does not imply that the two equations should be satisfied at the same time. They mean 
we can cut off the series in two different ways; while each cut off of the series provides a kind 
of solutions for the equation ( ^.26 ). We might have two kinds of solutions in this case. 



Case (2): In this case, let ^ < 1, so p = ^s(l - f ) > 0. From equations (Q, (|Aj| , ( |A?7| ) 
and ( [A.g| ), we can infer that Eq. ( [A.7] ) has no solutions for any t; while equations ( [A.5[ ), ( |A.6| ) 
and ( [A.8| ) might have solutions when t < {-^s^)"^. 
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In this case, f{t) is a decreasing function about t when t < (^sf ) ; while g{t) is still a 
increasing function about t when t < {^s^)'^. We obtain 



m / /aN 2 



By these equations, we determine the extent of n to be 
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For ^ < 1, 



+(l-f) 



< 0. So Eq. ( [A.20D is impossible, and the correspond- 



ing equation (A. 8) has no solutions. While Eq. ( [A. 181 ) is also impossible obviously, and the 



corresponding equation (A^) has no solutions. 

In summaries, in this case, only Eq. ( |A.19 ) is possible, and the corresponding equation 
( |A.6| ) has solutions. We only have a way to cut off the series, and we can have one kind of 
solutions corresponding to this cut off when f < 1. 

In addition, we make more discussions about a special case. We analyze whether equations 
( |A.5| ), ( [A.6| ), (A. 7) and ( [A.8| ) can have the solution t = 0. Replacing t with in these equations, 
we find that only Eq. ([A.5|) is possible. We obtain the condition 
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io V 6 



(A.21) 



For T > 1, this condition can be satisfied. When the condition Eq. ( A.21 ) is satisfied, the n' 



Jh 



massive solutions will coincide with the zero mode solutions given by Eq. (B.l) in appendix p[ 



B. Explicit solutions for zero modes and massive modes 



In this appendix, we give the solutions for zero modes and massive modes explicitly. In the 
following we discuss two cases: m > and m < 0. 

Case (1): For the positive bulk mass parameters, that is, m > 0. The normalizable zero 
mode for the metric (2.46) is given by 

Fo{x) = 0, Goix) 



UJ ma. .7 

X~~''b{x+ 1)" 



(B.l) 



where we have defined x = 
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For the massive modes, the solution determined by the parameter set ( 2.55D , ( |2.56| ) and 
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( |2.57D is given by 
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The solution for Gi{x) is determined by the equation 



Gn{x) = — 
An 



m X + X 
— s 



d 



, , fF„(x) -x-f Fn(a;) 
u; X + 1 ax 



n = 1,2,3, 



(B.2) 

(B.3) 
(B.4) 

(B.5) 
(B.6) 

(B.7) 



According to the discussions in subsection |2.3| , the massive modes emerge in pairs. The 
other solution in pairs with the solution (Fi, Gi), that is, the solution corresponding to the 
eigenvalue — Ai, is given by 



F^i(x) =Fi(x), G_i(x) = -Gi(x). 



(B.8) 



Case (2): Now we discuss the case m < 0. Redefine m = —m, then m > 0. The action 
becomes to be 



(B.9) 



After KK decompositions like in subsection 2.1, the equations ( p. 11 ) and (2.12) change to 
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Prom equations ( |B.13| ) and (|B.14 ), we see that Gn{z) conform to the similar equation Hke Fn{ 
in equations ( 2.26| ) and ( [2.271 ). The massive solutions for Gn{z) are given by 



G{z) = Cie"^'^^(e'^^ + 6)''-^hypergeom { p - ^ + 1^,1 - p - ^ + - 2^, 



+ C2e>^^{e'^'' + 6)"''"''hypergeom [ p + p + u,l - p + p + u;l + 2p, 



,(B.15) 



where p = f .(1 - f ), M = V(f (§)' " (^)'' = V(^s)' - Therefore, the 

conclusions about the extent of n in subsection ^ and appendix ^ applies here. We just need 
replace m with m, because we assume m > in those analysis. However, the solutions of zero 
mode and massive modes change. 

The zero mode solution is given by 



Fo{x) 



The massive mode for Gi{z) is given by 



Gi{x) 



Pi X 

71 X + 1 



in which 



7n 



l-2//„, n = l,2,3,' 



Pn = 1- p- Pn + l^n, p = —S [1 




bJ ' 



While the solution for Fi (x) is determined by the equation 

UJ 



Fn{x) 



An 



''-Gnix) + X — Gn{x) 

ax 



, n = l,2,3,' 



u X + 1 

For the other solution in pairs with the solution {Fi, Gi), now we obtain 

= -Fi{x), G_i(x) = Gi{x). 



(B.16) 

(B.17) 
(B.18) 

(B.19) 
(B.20) 

(B.21) 
(B.22) 

(B.23) 
(B.24) 



In order to determine the normalization constants Nq and A'^i in the above equations, we 
designate the normalization conditions as 



/ 



dz {F*Fn + GlGn) = 6nn, n = 0,1,-1, 



(B.25) 



where (5„n = 1 and n is not be summed. 



- 44 - 



C. Definitions for quantities 



We define the quantities in equations ( [4.1[ ) and ( |4.2[ ) as follows. For the field Q, we define that 

aQ = -il dziF^F^ + G^G^), bQ = j dz{F^*F\ + G^^), (C.l) 



Note that AA^ is pure imaginary according to our designation in subsection 4.1, so (Zq and ftj; 
are real. For the field P, we define that 



dv 



j dz{F^*Ff + G^*G^), bv = j dz{Ff*F'^^ + Gf (C.2) 



-aijAA^. 



Note that here AA^ is real according to our designation in subsection 4.1. Functions Fn{z) 
and Gn{z) are determined by equations (2.79) and ( p.SOp . They can be determined according 



to our discussions in subsection |2.5| . Obviously they have the similar forms and characters to 
the solutions given in appendix Due to the special characters of these solutions discussed 



in subsection |2.3| , the matrix K and other matrices in subsection all have very concise 
structures. 

For the quantities in Eq. (|4.3D , we have defined them as follows. For , we define that 



lb 
d 



dzB{z){F^*G'^ + GTF^), a 
dzB{z){F^*G'E + GfF^), c: 



dzB{z){F^*G'^^ + GfF^^). 



j dzBiz)iF^*GY + G^*F^), 
I dzB{z){F^*Gf + GfFr), 



For y^, we define that 



= I dzB{z){F^*Ff - G^*Gf ), h' = j dzB{z){F^*F^ - G^*G^), 
' = J dzB{z){F^*Ff - Gf), d' = y dzB{z){F^*F'E^ - G?*G?i^ 



The elements in equations ( |C.3| ) and ( p.4| ) can be complex numbers generally. 
The eigenvalues in equations (|4.4| ) and (^]^) are given by 



Af 



6q _ 1 
2 2 

bv_l 
2 2 



Sal + bl, A^ = l 



(C.3) 



(C.4) 



(C.5) 
(C.6) 
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The elements in equations ( [4. 61) and (4/7) are given as follows. For the field Q, we define 



that 



'T2 



)AX^ cos% + 6oAA^sin( 



bo = I -^QQAA^sin^ + 6qAA^cos( 



For the field P, we define that 

1 

.72 



ax> = [ —^ai)AX^ cos 9v — AA'^ sin 9x) 



1 

71 



aj}AX^ sin 9j} + AA^ cos 9j} 



V 



1 1 

V 1 V 3 

1 1 



1 1 
1 1 

/a?./a?' 



(C.7) 



(C.8) 



The elements of H'^'^ in Eq. ( ^.lOP are given as follows, sq, aq and Ib^ are given by 



: cos an cos '< 



: sm On sm 1 



(C.9) 



1/1 ^ „ 1 

sin 9qi cos vo H ^= cos t>o sm 1 



IbQ = — ^ I — ; — COS 9o sin i?© H sin cos -i 

v2 



While (Cq and dq are given by 
1 1 

<do = - 



: sin 9o sin i?o H ; — cos tfn cos 1 



1 1 

: sin sin t?o H ; — cos 6*0 cos 1 



(C.IO) 



The elements of H'^ in Eq. ( p^ ) are given as follows, s^i a© and bx) are given by 



1 1 

cos 9d cos i?!) ; sin 0x> sin I , 



(C.ll) 



ax) = —pz \ — ■ cos 9-1) sin t?© H ; sin 9^ cos 



hj) = ( — ■ sin 6'x) cos i?© H ; cos 9x) sin t^x" 



V2 \ /Xf 



While <cd and dx) are given by 



: sin 9x> sin t?x> 



: cos Oj) cos t?D 



(C.12) 



+ ( — ; sin 6x> sin i?x) 



: cos 6d cos i?D 



In order to make single- value decomposition for and Y^, we define two matrices as 
follows 



yDyD^ 



A* B C 
-A* C B 



( S' Al -Al 
B' C 
\-A!* C B' 



(C.13) 



The matrices M and M' are defined by the above matrix multiplication. Their elements can 
be expressed with the elements of and respectively. We omit their explicit expressions 
here. By the definitions in Eq. ( p.l3| ), the elements S, B, C, S' , B' and C are all real numbers. 
By the elements of these two new matrices, the quantities in equations ( [4.14 ) and (4.15) can 
be defined as follows. For y^, we define that 



{B-C + S) + \/(i3-C-cS)2 + 8|^|2 



(C.14) 



= - [(e - C + 5) - V(^-C-cS)2 + 8|^|2 



, s^ = e + c, 



^=1^1 exp iSy , tan ( 



'y 



{B-C-S) + y/{B - C - 5)2 + 8 I ^ |2 
2^/2 \A\ 



For Y^, we define that 
1 r 



1 



(B' -C' + S') + -C - S'f + 8 I ^' |2 , 

{B' -C + S') - V(S' - C' - cS')2 + 8 I A' |2j , = B' + C, 

(B' -C - S') + ^(^' - C - 5')2 + 8 I A! |2 



(C.15) 



= I I expi(5y, tan^y 



2^/2\ A! 
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